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Measurement of r polarization 
in e + e~ annihilation at \/s=58 GeV 
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The polarization of r leptons in the reaction e + e~ — ► r + r _ has been measured using a e + e~ 
collider, TRISTAN, at the center-of-mass energy of 58 GeV. From the kinematical distributions 
of daughter particles in r — > ez/P, fivu, pv or tt{K)u decays, the average polarization of r~ 
and its forward-backward asymmetry have been evaluated to be 0.012±0.058 and 0.029±0.057, 
respectively. 
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1 Introduction 



The polarization of r leptons in the reaction, 

e + e~ — > t + t~ (1) 

provides us with information concerning the properties of the neutral current of leptons [[jj . The 
measurement of it is, therefore, important in every Q 2 range of the reaction. 
The average polarization of r leptons is defined as 

= a(h = +l)-a(h = -l) 
/r a(h = +l) + a(h = -iy {) 

where h(= ±1) is the helicity of r~ . Another relevant observable is the forward-backward 
asymmetry of the polarization defined as 

u f Vf ~ <? b Vb 

Av = ^ j_„ ' ( 3 ) 

& F + &B 

where Vf and Vb are the average polarizations measured in the forward and backward regions, 
respectively, relative to the incoming e~ direction. 

Although there are many precise measurements on V T and A-p at Z° pole Q, very few are 
available below the pole (3|. Besides, the previous measurements below the pole suffered from 
relatively large statistical errors due to the small cross section of reaction (|l]). 

This letter describes a measurement of V T and A-p at a center-of-mass energy (\/~s) of 58 
GeV using the VENUS detector at TRISTAN. The analysis is based on the data corresponding 
to an integrated luminosity of 271 pb _1 . 

In the framework of the Standard Model [Q] , the differential cross section of the reaction (|l|) 
can be expressed for unpolarized electron and positron as 

^-(cos9,h) = ^ \F (s)(l + cos 2 6) + 2F 1 {s)cos8-h{F 2 (s){l + cos 2 6) + 2F 3 {s)cos6}} , (4) 
ail 8s L J 

where 9 is the scattering angle between e~ and t~. The four form- factors are given as 

F (s) = 1 + 2g Ve g VT Rex{s) + {g\ e + ff| e )(#y T + 9at)\x{s)\ 2 , 
Fi(s) = 2g Ae g AT Rex(s) + ^g Ve g Ae g VT g AT \x{s)\ 2 , 
F 2 (s) = 2g Ve g Ar Rex(s) + 2{g\ e + g Ae )gvTgAr\x(s)\ 2 , 

^30) = ^gAe9VT^ex(s) + 2g Ve g Ae (gy T + g 2 Ar )\x(s)\ 2 . (5) 

Here, gyi and g A \ are the vector and axial-vector coupling constants, respectively, of leptons I 
to Z°. The function x( s ) can be written using the mass of Z° (Mz), its width (T^) and the 



weak mixing angle (9w) as 

X ^ = 4 sin 2 9 W cos 2 9 W X s- M 2 z + isTz/Mz' ^ 
Integrating Eq. (|j) over the full solid angle, the polarization and its forward-backward asymme- 
try can be expressed as 

V = 

= "it" ^ 

Since the term including \x(s)\ are small at the TRISTAN energy, V T is sensitive to gvedAr, 
and to gAe9Vr- On the other hand, they are sensitive to A T and A e , respectively at the Z° 
pole, where A\ is defined as A\ = IgviQAil '{dvi + 9ai)- ^ n * n ^ s sense > the present measurement is 
complementary to those from the experiments at the pole. 



2 The VENUS Detector 

Since the VENUS detector is described in detail elsewhere ||], we briefly summarize those 
components relevant to the present measurement. The central drift chamber (CDC) is the main 
component for charged particle tracks, located in a uniform magnetic field of 0.75 T parallel to 
the beam axis produced by a superconducting solenoid. Charged particles with | cos#| < 0.75 
(9 is the angle with respect to the beam axis) are detected in the CDC with a momentum 
resolution of c p /p ~ 0.008pt(GeV/c), where pt is the transverse momentum with respect to 
the beam axis. The angular resolution is 8sin 2 # mrad and 1 mrad for polar and azimuthal 
angles, respectively. A large cylindrical transition radiation detector (TRD) surrounds the CDC 
and covers an angular region of | cos#| < 0.68. It provides an e/V discrimination capability as 
demonstrated in Fig. 1. 

Between the solenoid coil and the TRD, 96 time-of- flight (TOF) counters are located at a 
radius of 1.66 m. A time resolution of 200 ps is obtained for charged particles within | cos0| < 
0.78. A cylindrical array of 5160 lead-glass blocks (LG) is located outside the solenoid in order 
to measure the energy of electrons and photons within | cos 9\ < 0.80. Its good energy resolution 
(3.8% for 30GeV electrons) makes it possible to identify electrons and tt° clearly. Outside the 
iron return yoke, eight layers of muon chambers (MU) are stacked, of which six inner layers are 
aligned along the beam direction, while the outermost two layers are stacked orthogonal to the 
inner ones. In the present analysis, the inner six layers are used for muon identification in an 
angular range of | cos#| < 0.5. 
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3 Selection of r decays 

The selection of the r decay of a specific mode proceeded in two steps. First, an preselection 
was made to reject non-r events such as two-photon collision (77 events) and multihadron events. 
The r pair events are characterized with such topology as low-multiplicity, back-to-back narrow 
jets. The criteria of the preselection are as follows: 

1) The number of CDC tracks should be between two and eight, where only those tracks with 
Pt greater than 0.2 GeV/c coming from the interaction region are counted inside | cos6*| < 0.50. 

2) The total visible energy, the sum of the track momenta and the calorimeter energies must be 
greater than 0.2-^/5. 

3) The total energy deposited in the calorimeters must be smaller than 0.8^. 

4) The invariant mass calculated inside each of two hemispheres with respect to the thrust axis 
must be smaller than 4.0 GeV/c 2 . 

5) The acollinearity angle between the jet axes, defined by the vector sum of track momenta 
and energy flows in each hemisphere must be smaller than 40°. 

The conditions were effective to reject 77 events (2 and 5), Bhabha events (3) and multihadronic 
events (4). 

In the second step, decay modes into a single-prong final state were identified independently 
in each hemisphere. The decay modes considered here are r — > evu, puis, tt(K)u and pu, the 
sum of which amounts to 75 % of all r decays. 

r — ► pvv : The muons were identified by using the MU. The tracks which were associated 
with at least four corresponding hit layers of the MU were identified as p. Those events 
in which both hemispheres were identified as pvv were discarded to reject ee — ► pp and 
eepp events. 

t — ► pv '. For the non-^ tracks, the energy deposit in the TRD (£^trd) was required to be 
smaller than 10 keV to ensure the track to be ir^. Those tracks which were associated with 
one or two neutral energy clusters in the LG inside a 30° cone were taken as candidates 
of p. Figure 2a) shows the invariant mass distribution of the selected two clusters. A 
peak corresponding to 7r° is clearly seen above small background. We selected those pairs 
whose invariant mass lies within 35 MeV/c 2 around the 7r° mass as candidates. Any single 
isolated cluster was taken as ir° whose decay gammas merged into one cluster. Finally, 
the invariant mass of the 7r ± and the 7r° candidates were calculated. Its distribution is 
shown in Fig. 2b). We selected those which lie within 230 MeV/c 2 around the p mass as 
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pv decay candidates. 

r — ► evv : The electrons were identified by using the LG and TRD information ||. Among 
the tracks not selected in the preceding identifications, the candidate tracks were selected 
by requiring E/p, the ratio of energy deposit in the LG to the track momentum, to be 
greater than 0.8. -Etrd was also required to exceed 10 keV to assure the passage of 
electron. Those events in which both hemispheres were identified as a evv decay were 
removed in order to reject Bhabha events and ee — > eeee events. 

r — ► tt{K)v : Among those tracks which were associated with no substantial MU hits (one hit 
layer was allowed), tracks were selected by requiring E/p < 0.8 and £trd <10 keV. 
They were also required to have no neutral clusters with energy greater than 200 MeV 
within 60° around them. Since tt/K separation was not possible in the present analysis, 
r — > Kv mode was also included in this category. 

Application of the TRD, which is one of the unique features of the VENUS detector, assures 
reliable identification among evv, pvv and ttv decay modes. 

The number of the selected decays and the estimated background are summarized in Table 
1. The efficiencies and contaminations for the above identification procedures were evaluated as 
a function of the daughter particle momentum, using real data from other processes: ee — > ee, 
ee7 and eeee as the electron sample; ee — > p,p, fip'y, eep,p and cosmic rays as the muon sample; 
ee — > eeim, r — > 3ttv and K s — * 2ir as the pion sample. It was found that there were no apparent 
momentum dependence in the efficiencies nor contaminations. They are also summarized in 
Table 1. For a confirmation, the branching ratios evaluated from our data are also listed together 
with the world average Those are in good agreements with one another. 



4 Evaluation of the polarization and its asymmetry 

In the case of the three body decays of r leptons such as r — * evv or fivv, the energy 
spectrum of the charged leptons is expressed as || 

^- = a {x)+V T b{x), (8) 

where a(x) = (5 - 9x 2 + Ax 3 )/3, b(x) = (1 - 9x 2 + 8x 3 )/3 and x is the laboratory energy of 
the charged lepton scaled to the beam energy. Figure 3 shows the observed x distributions for 
&)evv and b)pvv modes. 
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For the two body decays such as r — > ttu or pv, the momentum spectrum of the decay 
particle is expressed in the same way as || 

= i + aVr {2 x - l) (9) 

where a is an analyzing power related to the spin of the final state (a=l for ttv mode and a=0.46 
for pv mode). Figure 3c) shows the observed spectra for ttv mode. The smaller analyzing power 
for the pv mode, which comes from the mixture of opposite polarization states of p, can be 
restored by introducing a second kinematical variable || to discriminate the polarization state. 
In the present analysis, we used the decay angle of p in the r rest frame relative to the r flight 
direction (cosip T ~ 2x — 1) and the decay angle of the charged ir in the p rest frame relative 
to the p flight direction (cosip p ). The cosf/'p distribution was examined in two cos^V regions, 
cost/v < and cos^V > 0, separately. The observed distributions are shown in Fig. 4. 

The polarization was evaluated by fitting the linear combination of the kinematical dis- 
tributions expected for the full polarization cases {h = ±1) to the corresponding observed 
distributions. The fit was done by minimizing a x 2 defined as 



x 2 = E 



jymeas _![(! + V T ){N+ + Nf G+ ) + (1 - V T ){Nr + jVf G " 



(10) 



where N™ eas is the number of entries in i-th bin, N^~(N~) and N® G+ (N® G ~) are those of 
Monte Carlo events with h = +1(— 1) for the signal and the background coming from other 
decay modes of r, respectively. The denominator, Oi is the statistical error in each bin. The 
expected contributions from non-r process were subtracted from NF 1 ^ beforehand. 

The Monte Carlo events for each helicity state of r leptons were generated by using KORALZ 
4.0 [10] for r pair production and TAUOLA2.5 [ll]] for r decay. The radiative effect in the initial 
state, which affects the energy spectra substantially, is properly incorporated up to 0(a 2 ) in the 
event generation. The efficiencies in the decay mode identifications were treated according to 
the functions evaluated in the previous section. The results of the fit are shown in Figs. 3 and 
4 with histograms. The average polarizations of r thus evaluated are tabulated in Table 2. 

The forward-backward asymmetry of the polarization, A-p, was evaluated by the same pro- 
cedure as above in the forward (0 < cos 9 < 0.5) and backward (—0.5 < cos 9 < 0) region, 
separately. The scattering angle was determined by the thrust axis of the event and the charge 
of the daughter particle. From the forward and backward polarization thus evaluated, its asym- 
metry was calculated as listed in Table 3. 

The systematic errors in the above measurements were considered in three categories. 
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Detection efficiency: The evaluation of polarization can be affected by uncertainties in the mo- 
mentum dependence of the efficiency of the decay mode identification. Since the estimation was 
made using real data, there are substantial statistical uncertainties. Their effects on V T were 
estimated by changing the efficiency independently in each momentum bin within their errors. 
Background: Uncertainties related to the background, either from mis-identified r decays or 
non-r processes, can lead to the systematic error. The evaluation were found to be very stable 
with the possible variation of the background fraction. 

Monte Carlo statistics : Because of limited statistics in the Monte Carlo simulation, the gen- 
erated spectra had some uncertainties. The systematic errors due to them were estimated by 
changing the number of entries in each bin of the spectra within their statistical errors. 
The systematic errors thus estimated are summarized in Table 2 and Table 3. As found in the 
tables, the systematic errors are significantly smaller than the statistical ones in the present 
analysis. 

Combining the results for all the modes studied in the present analysis with neglecting the 
small systematic errors safely, the average polarization of r leptons and its forward-backward 
asymmetry are 

V T = 0.012 ± 0.058, (11) 
A P = 0.029 ± 0.057. (12) 

Although the present results are consistent with zero and are also consistent with the values 
predicted by the Standard Model (0.028 and 0.021, respectively, at y/s = 58 GeV), they are 
the first substantial measurements in the energy region below the Z° pole as found in Fig. 5, 
where the results from the present and the previous experiments are plotted together with the 
standard model prediction. 

5 Conclusions 

The polarization of r leptons and its forward-backward asymmetry in the reaction e + e~ — ► r + r 
have been measured at the center-of-mass energy of 58 GeV. They are evaluated to be 0.012 ± 
0.058 and 0.029 ±0.057, respectively, combining the results for all the decay modes studied. The 
prediction of the standard model is confirmed in the energy region below the Z° pole. 
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Table 1: The number of identified decays, the estimated background contribution and the 
selection efficiency for each decay mode. Also listed are calculated branching ratios and 
corresponding PDG96 values. 





pv 


n(K)v 


evv 


pVV 


# observed 


829 


287 


564 


628 


# background(r) 


75.6 ±8.1 


25.9 ± 1.6 


2.0 ±0.4 


6.2 ±0.8 


# background(non-r) 


17.7 ±2.4 


8.7 ± 1.6 


20.3 ±4.0 


21.7 ± 1.9 


Efficiency(%) 


36.1 ± 1.0 


24.3 ±0.8 


34.5 ± 1.0 


41.3 ± 1.3 


BR(%) 


24.1 ±0.9 ±0.8 


12.3 ±0.8 ±0.5 


18.6 ±0.8 ±0.7 


17.2 ±0.7 ±0.6 


PDG96(%) 


25.24 ±0.16 


12.03 ±0.14 


17.83 ± 0.08 


17.35 ±0.10 



Table 2: Evaluated V T and systematic errors contributing to it. Errors indicated in the 
first column represent statistical ones. 





evv 


pvv 


n(K)v 


pv 




0.13 ±0.18 


0.15 ±0.20 


-0.01 ±0.12 


-0.03 ±0.07 


Efficiency 


0.021 


0.021 


0.018 


0.024 


r BG 


0.014 


0.013 


0.010 


0.009 


Non-r BG 


0.012 


0.008 


0.009 


0.011 


Monte Carlo statistics. 


0.021 


0.020 


0.013 


0.015 


Total systematic error 


0.035 


0.033 


0.026 


0.032 
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Table 3: Evaluated A-p and systematic errors contributing to it. Errors indicated in 
first column represent statistical ones. 





evv 


\xvv 


ir{K)v 


pv 


A P 


-0.02 ±0.17 


-0.01 ±0.19 


0.03 ±0.12 


0.04 ± 0.07 


Efficiency 


0.023 


0.025 


0.017 


0.026 


r BG 


0.015 


0.012 


0.011 


0.009 


Non-r BG 


0.013 


0.008 


0.010 


0.014 


Monte Carlo statistics. 


0.022 


0.019 


0.015 


0.013 


Total systematic error 


0.038 


0.035 


0.027 


0.033 
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Figure 1: Energy deposit spectra in the TRD, for electrons (solid line histogram) and for 
muons (hatched histogram) with their momenta greater than 2 GeV/c. The spectrum 
for pions must be the same as that for muons. The cut at 10 keV provides a clear 
discrimination between electrons and pions. 
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Figure 2: Invariant mass distribution of a)the two photons and b)^^ in the selection 
of the pv decay mode. The solid histogram shows the simulation result. The estimated 
background contributions are indicated with the dashed histogram. 
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Figure 3: Scaled energy distribution for a)r — > ez/z/, b)fivv and c)7rzv modes. The filled 
circles indicate the selected candidates. The solid line histograms represent the result 
of the fit. Contributions from each helicity state are indicated by the dashed (h = +1) 
and dotted (h = —1) histograms. The background contributions are illustrated with the 
shaded histogram. 
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Figure 4: Observed cos ip p distributions for the r — > pv candidates in the region of a)0 < 
cosi\) T and b)cos-0 T < 0. The filled circles indicate the selected candidate. The solid 
line histograms represent the result of the fit. Contributions from each helicity state are 
indicated by the dashed (h = +1) and dotted (h = —1) histograms. The background 
contributions are illustrated with the shaded histogram. 
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Figure 5: Measured a)average polarization V T and repolarization asymmetry A-p together 
with those from other experiments. Solid line indicates the Standard Model prediction. 
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